Distance matrices and isometric embeddings 



E. Bogomolny, O. Bohigas, and C. Schmit 
Universite Patis-Sud, CNRS, UMR 8626, 
Laboratoire de Physique Theorique et Modeles Statistiques, 
91405 Orsay Cedex, France 

February 2, 2008 
Abstract 

We review the relations between distance matrices and isometric embed- 
dings and give simple proofs that distance matrices defined on euclidean 
and spherical spaces have all eigenvalues except one non-positive. Several 
generalizations are discussed. 

1 Introduction 

Matrices with random (or pseudo-random) elements appear naturally in different 
physical problems and their statistical properties have been thoroughly investi- 
gated (see e.g. Ilj). A special case of random matrices, called distance matrices, 
has been recently proposed in [2]. They are defined for any metric space X with 
a probability measure /x on it as follows. Choose points Xj E X randomly 
distributed according to the measure fi. The matrix element M,j of the N x N 
distance matrix M equals the distance on X between points Xi and Xj 

Mij = distance(xj, Xj), i, j = 1, . . . , A . (1) 

In [3] we have discussed the eigenvalue density for distance matrices defined on 
certain manifolds. When first numerical calculations were performed, an intrigu- 
ing fact was observed, namely, all eigenvalues (except one) of distance matrices on 
euclidean and spherical manifolds were non-positive. However, this property was 
not fulfilled e.g. for points on a torus. Typically, eigenvalues of generic random 
matrices occupy the whole available energy space and to impose the condition 
that all of them but one are non-positive requires the control of the signs of all 
principal minors (see Section [S]) which is usually difficult to impose. In inves- 
tigating this fact we have found a direct proof that for distance matrices over 
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manifolds embedded into the euchdean space this property is automatically ful- 
filled. This relation between very basic geometrical properties of a manifold and 
spectral properties of its distance matrix was unexpected for us but analysing 
the literature we found that it has been proved by Schoenberg in the thirties 
[U \5\ and in [3] we noted this fact without details. Nevertheless, after many 
discussions on different occasions it became clear that this type of problems is 
practically unknown in the physical community and we think that it is of interest 
to present simple proofs of the main statements. The material in this note is not 
new (general references are [6]- [8]) but it seems that it has not been discussed in 
the random matrix community. 

By definition of distance the matrix elements of a distance matrix have the 
following properties 

a) positivity 

Mij > and Mij = only when i = j , (2) 

b) symmetry 

= Mji , (3) 

c) triangular inequalities 

Mij < Mik + Mkj for all i, j, k . (4) 

Eigenvalues, Ap, and eigenvectors, u'f^, of distance matrices are defined in the 
usual way 

N 

Y: M,,uf = K,v!f^ . (5) 
i=i 

Distance matrices ([TJ are real symmetric matrices and their eigenvalues are real. 
As all matrix elements of distance matrices are non-negative, the application of 
the Perron- Frobenius theorem ([9j V. 2 p. 49) states that these matrices have one 
special positive eigenvalue Aq > with the largest modulus. All other eigenvalues 
obey the inequality 

|Aj|<Ao. (6) 

As distance matrices have only real eigenvalues the equality is possible only if 
there is a negative eigenvalue A' = — Aq. 

The subject of this note is to demonstrate that eigenvalues of distance matrices 
defined on the euclidean or a spherical spac^ are all non-positive 

Ai<0, i = l,...,N-l, (7) 

^It means that the points xj lie in the d-dimensional euchdean space or on a sphere. 
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except the above-mentioned Perron-Frobenius eigenvalue Aq and that this re- 
markable property is mainly a consequence of the possibility of isometric em- 
bedding of a finite metric space with a given distance matrix into the euclidean 
space. 

We also remark that if, instead of the distance matrix ([T]) , one considers new 
matrices 

M^J^ = [distance(xi,fj)]'^ , (8) 

their eigenvalues also obey inequality ^ provided the exponent in the range 
< 7 < 2 for the euclidean space and < 7 < 1 for the spherical one. 

The plan of this paper is the following. In Section [2] it is demonstrated that if 
a finite metric space can be isometrically embedded into the euclidean space, then 
the matrix whose elements are the squares of distances between initial points is of 
negative type (cf. (j20p - (j2ip ). The inverse theorem is also true, namely, if a matrix 

is of negative type then the matrix with elements yJNij can be embedded into 
the euclidean space. A direct proof of the main theorem that matrices of negative 
type have all eigenvalues, except one, non-positive is presented in Section [3J In 
Section Slit is demonstrated that if a matrix Nij is of negative type, a new matrix 
with < 7 < 1 will also be of negative type. The general form of such metric 
transforms is also shortly discussed in this Section. In Section [5] spherical spaces 
are discussed and in Section [6] a simple proof that geodesic distance matrices for 
the spherical spaces are of negative type is presented. A resume of the results is 
given in Section [71 The derivation of the Cayley-Menger formula for the volume 
of a multi-dimensional simplex is reproduced for completeness in the Appendix. 



2 Isometric embedding 

Assume that we know a finite matrix M whose matrix elements M,-,- = 
1,...,A^) obey all properties of a distance (HJ-dH). The isometric embedding 
into the euclidean space consists in finding points Xi, if any, belonging to an 
euclidean space BP" such that the euclidean distance between each pair of points 
i,j coincides with Mij 

\\xi - XjW = Mij , (9) 
for all i, J = 1, . . . , A^. Here 1 1 . . . 1 1 is the euclidean distance 



\ fe=i 



(A;) 

and x\ with k = 1,... ,n are the euclidean coordinates of the n-dimensional 



point Xi 
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The necessary and sufficient conditions of the existence of solutions of Eq. ([9]) 
can be obtained from the following considerations [1] and [6l[7]. Choose a point, 
say xn and consider the vectors yi = Xi — xn with i = 1, . . . , N — 1. They form 
a simplex in the n-dimensional space i?". Construct the (A^ — 1) x n matrix of 
coordinates of these vectors 

V^k = yl''^ i = l,...,N-l, k=l,...,n (11) 

and multiply it by its transpose. The result is a (A^ — 1) x (N — 1) real symmetric 
matrix C = V'^ ■ V of scalar products 

Cij = yi-yj, i,j = l,...,N -1 . (12) 

Because vectors yj belong to the euclidean space their scalar products can be 
expressed through the distances between points 

Vi ■ Vj = ^(llyilP + llyjiP - Wvi - ^ilP) • (i3) 

Therefore the matrix Cij can be calculated from the squares of matrix elements 
of the distance matrix 

Q, = ^(M2^ + M|^-M2.) . (14) 

If points Xj obeying Eqs. ([9]) do exist then by construction the matrix Cij is such 
that the quadratic form 

7V-1 / N 

(15) 



N-l / N 

i,i=l \j=l 



is non-negative (> 0) for any choice of real numbers ^i,^2) • • • ,(,N-i- Inversely, if 
one has a symmetric positive matrix C, it can be written in the form 

C = U^U, (16) 

where the matrix U can be chosen, e.g., in the lower triangular form (the Cholesky 

(k) (k) 

decomposition). Then the elements of C/ = y!- give directly the coordinates 
of points obeying (jl4p which solve the problem of the isometric embedding. 

The quadratic form (jl5p can be rewritten in a simpler form by introducing a 
new variable £,n = — J2f=i^ ■ Then 

1 ^ 

(^^0 = -2 E ^U^^i ■ (17) 
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Therefore the necessary and sufficient condition of the existence of isometric 
embedding of a finite metric space with the distance matrix M into the euchdean 
space is that a new matrix N whose matrix elements equal the square of matrix 
elements of the matrix M 

N,, = M2. (18) 
is such that the quadratic form associated with it 

N 

(W) = E ^iMj (19) 

is non-positive 

TV 

for all choices of real numbers S^j, j = 1, . . . , N with zero sum 

N 

i=i 

In general, a real symmetric matrix obeying these conditions is called a matrix 
of negative type. 

The Schoenberg theorem states that if a metric space with a distance matrix 
Mij can be isometrically embedded into the euclidean space, the matrix Mfj is of 
negative type and if a matrix Nij is of negative type, the metric space with the 
distance matrix ^JNij can be isometrically embedded into the euclidean space. 
The minimal dimension of the embedded euclidean space is the rank of the matrix 
Cij in Eq. ([HI). 

3 Eigenvalues of negative type matrices 

In this Section we present, following [5j, the direct proof that any matrix Nij of 
the negative type has all eigenvalues except one non-positive (< 0). An indirect 
proof of this statement can be found in [8] . 

The law of inertia (see e.g. [9] V. 1 p. 298) states that if a real quadratic form 

N 

{xNx) = E NijXiXj (22) 
is transformed into a sum of squares of independent linear forms Xi 

r 

(xNx) = E b^Xf 

i=l 



2.7 = 1 ' 



(23) 
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then the total number of positive and negative coefficients 6j is independent of 
the representation. In particular, in the eigenbasis of the real symmetric matrix 

N, 



N 



(xNx) = J2^i 



(24) 



1=1 



and the law of inertia permits to determine the number of positive and negative 
eigenvalues Aj. 

According to the Jacobi theorem (see e.g. p_j V. 1 p. 305) if the principal 
minors Aj of a matrix are non-zero then the number of positive (resp. negative) 
terms in (I23p coincides with the number of conservation (resp. alteration) of signs 
in the sequence 

l,Ai,A2,...,A;v (25) 

(we assume that the matrix a,-,- is of full rank). Recall that the principal minor 
An of a matrix N is the determinant of the left-upper n x n sub-matrix 



An 



Nn Ni2 iVi3 
Ni2 N22 N23 

Nm N2n ... 



n— 1 n 



N2n 
Nnn 



(26) 



For distance matrices Nn = and Ai = which prevents the direct application 
of the Jacobi theorem. This formal difficulty can, for example, be overcome as 
follows. It is clear that eigenvalues and other principal minors of generic distance 
matrices are non-zero. Therefore if one adds to a diagonal matrix e5ij with 
small e the signs of eigenvalues will not change. But in such a case Ai = e and 
the sequence takes the form 



l,e,A2,A3 



.A 



N 



(27) 



We shall prove below that principal minors A„ of distance matrices of the negative 
type have alternating sign 



n 



2,3,... ,N . 



(28) 



Irrespective of the sign of e there is one conservation of sign and A^ — 1 alterations 
of signs in the sequence (j27p . Therefore, according to the Jacobi theorem, distance 
matrices of the negative type have one positive (the Perron-Frobenius) eigenvalue 
and all other eigenvalues are non-positive. 

Because the matrix A^j-,- is of the negative type, the metric space with the 
distance ^/Njl can be embedded into the euclidean space. It means that there 
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exist points xj in the euclidean space such that the euchdean distances between 
any pairs of points equal 



Di 



(29) 



Let us consider one of these points as the origin (say xi). Points X2, ■ ■ ■ ,Xn can 
be viewed as vertices of a (n — 1) -dimensional simplex. Denote Du = ri. Then 
the distance between any pair of points can be expressed as follows 



R 



+ r? 



2rjrj cos (pij 



(30) 



where ^pij is the euclidean angle between vectors Xi — xi and xj — xi. 

Let us perform an inversion — > 1/rj for all i = 2, . . . ,n. Then instead of 
n—1 points X2, • • • , we get a new set of n — 1 euclidean points X2, . . . , x„ whose 
mutual distances Dij can be expressed through the old distances as 



1 



1 



'12+12-2 



1 



nr. 



COS ipij 



DuD 



(31) 



Because the points Xj belong to the euclidean space the new points Xj with 
j = 2, . . . , N plus the point xi form an — 1-dimensional euclidean simplex. The 
volume of this simplex can be computed by the Cayley-Menger determinantal for- 
mula (see e.g. |12j p. 124 and also [T3] for an early reference) which expresses the 
volume V{Pi, . . . ,Pn) of a n-dimensional euclidean simplex through the lengths 
of its sides 



p 



-1 



2"-i[(n 

where the Cayley-Menger determinant is 



1)!]2 



(32) 



DiPi,...,Pn) 



1 





D2 
^12 



1 Dl, Dl^ 



1 





(33) 



For com- 



and Dij are the distances between points i and j for i^j = 1, . . . ,n. 
pleteness we give in Appendix a derivation of this formula. 

In our case the lengths of the transformed simplex are given by Eq. ()3ip . 
As each Dij = ^JNij , the squares of the lengths which enter the Cayley-Menger 
formula (l33l) are 



Dt 



NuNu 



(34) 
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Therefore for each n = 2, 



= 2,... 


N 











1 


1 




1 
1 




N32 


N23 
N12N13 




■Vn) = 


N13N12 




1 


N„2 






Ni„Ni2 





1 

N2„ 



Ni3Ni„ 





(35) 



As the determinant is a multi-hnear form of row and columns by multiphcation 
of each row (ij) and each column (ji) by N^j one gets 



Dim 



[Ni2N-_ 



13 



1N\ 



N12 Ni3 

N21 iV23 

N31 Ns2 



n2 



Nn3 



N2n 





(36) 



But the determinant in this expression coincides with the principal minor of the 
initial distance matrix. Therefore 

A„ = (-l)-"i[iVi2iVi3 . . . N,N?2^-'[{n - ly.fvHvi, • • • , ^n) , (37) 

which proves that the principal minors of matrices of the negative type are of 
alternate signs. This relation, as explained above, implies that all eigenvalues of 
such matrices (except possibly one) are non-positive. 



4 Metric transform 

The problem of isometric embedding gives rise to different generalizations. One 
type of question is the following. Let the points xj with j = 1, . . . , be points of 
the euclidean space i?". Find all functions F{r) (called metric transforms) such 
that the finite metric space with the distance matrix 

Mij = Fi\\xi - XjW) (38) 

can be embedded into an euclidean space with certain k. Here || 
the euclidean distance (fTOl) between point Xi and Xj. 

In [5] it was proved that general metric transforms can be expressed through 
radial positive definite functions. A real function /(r) is called radial positive 
definite provided 

N 
«J=1 
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for all choices of points Xj G i?" and of real numbers ^. 
An important example of such a function is 

/(r) = exp(-AV2) . (40) 

The positive definite property of this function is the direct consequence of the 
well known formula 

exp(-A2p||2) = ^-i-^ |^^e^^-^^exp(-P|p)d-fc , (41) 

from which it follows that 

N , „ 

-(4vr)"/2^ 

The following theorem is easily proved [5] . The finite metric space with a distance 
matrix Mij can be isometrically embedded into the euclidean space if and only if 
the quadratic form 

N 

expi-X'Mf^)Uj (43) 

is non- negative (> 0) for all choices of real numbers and all A — > 0. 

The proof is as follows. If the space can be isometrically embedded into the 
euclidean space then there exist points Xj E i?" such that Mij = \\xi — Xj\\. 
Because exp (-AV) is a radial positive definite function the quadratic form (|43p 
is non-negative. Conversely, if the quadratic form is non-negative for A — > then 

N 

5](l-AM2)e.e, >0 (44) 

i.j=l 

for all ^j. Choosing J^jLi^j = cancels the first term and reduces the above 
inequality to (jl9p . thus proving the existence of the embedding. 

The fact that exp(— A^r^) is a radially positive definite function permits also 
to prove [5] that the metric space with the distance equal a power of the euclidean 
distance 

Mlj = \\xi - Xj\r i,j = l,...,N, (45) 

where Xi G i?" and < 7 < 1 can be embedded into the euclidean space. The 
proof follows from the identity valid for < 7 < 1 

POO ^ \ 

lf" = c,y^ (l-exp(-A2t2))_^ ^ (4g) 



iXxi-k 



-\\k\? 



(i"A: > . 



(42) 
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with 

dX 



(l-exp(A2))^>0. (47) 



One has 

Ell-, ~* 1 1 2'7 c c V ^ ^ J. ^2 11^^ ^ . 1 1 2 , d\ 

\\Xi — XjW qic^j = Cj y ^ qic,j j [i — e )- 



fiiT.^^)' - (E«-'^"''-'^"'e.e,)]^ . (48) 



Al+27 

dX 

\ 2 2 

If = the first term is zero and as e~ ^ is radial positive definite, the 

right-hand side is negative which proves that the matrix 

I ^X'i j 1 1 



with < 7 < 1 is of negative type and the metric space with the distance ([451) 
can be embedded into the euchdean space. 

Combining together the above theorems, one concludes that if a matrix Nij 
is of negative type, then the matrix N^- with < 7 < 1 is also of negative type 
and all its eigenvalues, except at most one, are non-positive. 

General radial positive definite functions /(r) have the form [5] 

f{r) = / VLN{ru)d^ji{u) , (50) 

where the measure /i is non-negative, ii{u) > 0, and the function ^^{r) is the 
integral of e^^'^ with ||x|| = r over the (A^ — l)-dimensional sphere 

Here oj^-i = 2tt^^'^ /r{N/2) is the volume of the (A^ — l)-dimensional sphere, 
r(a;) is the Gamma function, and Jn{x) is the Bessel function. 

From the theorem ()43p it follows that the general form of a metric transform 

is 

^M = {f i^^<iK«)}"' (52) 

with a positive measure v{u) such that dv{u) /v?' exists. 
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5 Spherical spaces 



Eq. (j52p gives the general form of the metric transforms which transform an 
euchdean space into another euchdean space. Similar questions can be asked 
about the unit radius spherical spaced Sd-i which consist of points xj G 
obeying 

xl + xl + ... + xl = l . (53) 
The geodesic distance on the sphere is 

d{x, y ) = arccos(x • y ) ■ (54) 

The necessary and sufficient conditions that a metric space with the distance 
matrix Mij can be embedded isometrically into the spherical space with the 
distance ()54p coincide with the condition that initial points plus one point at 
the origin can be embedded into the euchdean space. From ()13p it follows that 
the later can be expressed as the non-negativity condition of the quadratic form 

N 

^ cos(M,,)e^O > (55) 

for all choices of real numbers ^j. 

Similarly as for the euchdean spaces one can find all positive definite functions 
on the spherical spaces. In [1^ it was proved that these functions have the form 

oo 

g{t) = J2 aiCf \cos t) , (56) 

1=0 

where all coefficients a/ are non-negative ai > 0. Here p = d — 2 and C/^(cost) 
are the Gegenbauer polynomials. 

This condition can easily be understood from the expression of the Gegen- 
bauer polynomial through the orthogonal set of the hyper-spherical harmonics 
y}''\x) (see e.g. lllj, 11.4.2) 

where h{p, I) is the dimension of the irreducible representations of the d — 1 
dimensional rotation group 

h{p,l) = {2l+p)^-^^^-^ . (58) 
p. t. 



^Modifications for splierical spaces of radius R are evident 
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If Eq. ([56]) is fulfilled, one has 

" ' N 



N oo h{p,l) ^ 



J2 9{d{x^,Xj))C^^j = 0Jd-iY.T7r-n 53 
i,j=i 1=0 "-^P' ^> k=l 

which is evidently non-negative (> 0). 



(59) 



6 Embedding of the spherical space into the euchdean 
space 

In this Section we show that distance matrices resulting from spherical geodesic 
distances are of negative type and, consequently, the metric space with the dis- 
tance equal the square root of spherical distances can be embedded into the 
euclidean space. 

The proof is based on a following lemma: the spherical geodesic distance (j54p 
has the expansion 

d{x, y) = arccos(f • y ) = Aq + X! ^«C'f^^(^ • V) > (60) 

«=odd 

where all with odd / are negative but Aq is positive. 

Eq. (|6Up is the expansion of arccos(i) over d-dimensional spherical harmonics. 
The coefficients A/ of this series are 

A, = r eCf^^ (cos 9) sinP ddO , (61) 

where hi (p) is the normalization integral of the Gegenbauer polynomials 

hi{p)= r [Cf {cos e)fsmPede (62) 
Jo 

whose explicit expression is (see e.g. |TT], 10.9.7) 



_ ^(;+p-i)!r((p + i)/2) 
- (/+p/2)/!(p-i)!r(p/2) • 



As Cn = 1, one gets 



Ao = I . (64) 

To compute Xi with / / it is convenient to use the Gegenbauer integral (see e.g. 
[n], 10.9.38) 

n\ j\'''°''C^{cose)s\n^''dde = 2^V^r(A + l/2) ^^[!^^^^^^ i"z-^Jn+A(^) , 

(65) 
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from which one obtains (cf. [TT], 11.4) 

Xi = i'2Pl\l+p/2)T{p/2) / t-Pl^Ji+^/2{t)f{t)dt , (66) 



where f{t) is the Fourier transform of the initial function 

fii) = ^ r esinee-^'^-^'de = - Jo{t)) . (67) 

Corresponding to the two terms in f{t) there are two terms in A;. The integral 
including e** is zero for all / 7^ and the integral with Jo{t) is zero for even /. For 
odd / 

Xi = -i''HPl\l+p/2)T{p/2) / r^-^'^Ji+j,/2{t)Mt)dt . (68) 







The last integral can be computed using the integral (flTj, 7.7.4.30) 



00 







t~PJ^{t)Mt)dt = (69) 

r(;9)r((i + + ^ - p)/2) 

2Pr((l + + p)/2)T{{l + 1/ + /i + ,9)/2)r((l + + p)/2) ' 

The final result is 



A; 



p{p + 21) 
Svr 



r(p/2)r(//2) 



r(i + (/+p)/2) 



2 



(70) 



This expression is negative which proves the lemma. 
Using this lemma and Eq. ([57|) , one concludes that 

N N . h{p,l) - 2 

E dix.,x,M, = Xoij:(.?+ E E 

i,j=l i=l /=odd ^' k=l 



N 



(71) 



As all A; with / > 1 are negative, this expression is negative for all choices of 



such that Y^f=i ij = 0, i.e. the spherical geodesic distance matrices are of 



negative type. 

From the theorem of the preceding Sections it follows that a new metric space 
with the distance 

d^'^\x, y) = [arccos(x • y )]'*' (72) 
is also of negative type when < 7 < 1 and the space with the distance 

[arccos(2; • y )]'''''^ (73) 

can be isometrically embedded into the euclidean space. 
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7 Conclusion 



The distance matrices for points in the euchdean and spherical spaces are of neg- 
ative type and, consequently, they have all eigenvalues, except one, non-positive. 

More generally, if points Xj belong to the euclidean space, the above statement 
is true for the matrices 

pi-fjipT (74) 

with < 7 < 1. 

If points Xj belong to the spherical space with the distance d{xi,Xj) given by 
Eq. (|54|) then the matrix 

(P{xi,Xj) (75) 

with < 7 < 1 is of negative type and has all eigenvalues, except one, non- 
negative. 

The following theorems are also of interest. 
The matrices with elements 

exp(— A^l |xj — Xjl P'^) (76) 

with < 7 < 1 are positive definite and have all eigenvalues positive for all 
A — > 0. For 7 = 1 this fact has been mentioned in |14j . 

The similar theorem for the spherical space states that matrices 

exp(— A^(i'''(xj, Xj)) (77) 

with < 7 < 1 are positive definite for all A. 
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Appendix 

The purpose of this Appendix is to give, following [12j . a proof of the Cayley- 
Menger determinantal formula ()32p . 

The volume Vn of the n-dimensional Euclidean simplex with one vertex on a 
point Xn+i and n vertices on points Xj with j = 1, . . . ,n is proportional to the 
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determinant of components of the n vectors xj — Xn+i 



nl 



„{i) 



X. 



(1) 

n+l 



X 



(2) 



(1) ^(1) 



„(2) 



(2) ^(2) 



n+l 



3/0 CC, 



n+l 



X 



X. 



in) _ in) 
1 -^n+l 
{n) _ (n) 
•^n+1 



„(1) ^(1) 



X, 



n+l 



(2) ^(2) 
■^n+1 



(n) 
Xn 



(n) 
''n+l 



(78) 



As above the subscripts denote the points and the superscripts denote their co- 
ordinates. This expression can be rewritten in a more symmetric form through 
the determinant of the (n + 1) x (n + 1) matrix 



Vn = -, 



X. 



(1) 



n+l ^n+1 



(2) 



.(2) 

^n 

(2) 



X 



X. 



(n) 
1 

(n) 



(n) 

(n) 
'^n+1 



(79) 



Simple manipulations show that it can be transformed in two different ways 



V . det An 



- — det B„ 
nl 



(80) 



where the (n + 2) x (n + 2) matrices An and Bn have the following forms 



/ 



An — 



1 

(xi)2 





-2x 
-2x. 



(1) 
1 

(1) 





-2x 



(n) 
1 
in) 



and 



B„ 



{Xn? 
V (^n+l)^ 

/ 

1 

1 , 



-2x1^) 



-2x 



-2x1") 



-2x 



(1) 

n+l 



-2x, 



(n) 
n+l 



1 

V 1 



X, 



-■1 



^n 

(1) 
n+l 



X 



X. 



(n) 
1 

in) 



An) 

„(«) 
"n+l 



1 

(^l)^ 



(Xn)2 
(^n+l)^ ) 



\ 

1 
1 

1 

\ 



(81) 



(82) 



Notice the position of the column of 1 in Bn- Therefore 

(-1)"+! 



' ft 



2^{n\f 



det Cn 



(83) 
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where C„ = An ■ B^. 

Direct calculations give the Cayler-Menger formula 



\n+l 



2'^(n!) 



1^2 



/ 


1 


1 


1 


1 







-'^1 n+1 


1 







-^2 n+1 


1 


^1 n 





-^n n+1 


V 1 


^1 n+1 


• • • n+1 






(84) 



where D 



3y 1 3^ -J 



is the length of the edge of the n-dimensional simplex. 
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